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Abstract: 

The amplitude of subthreshold meson production is calculated using dom- 
inant tree-level diagrams for three-body collisions. It is shown that the pro- 
duction can overwhelmingly be described by two-step processes. The effect 
of the genuine three-body contribution (i.e. the contribution which cannot 
be factorized) is discussed. The production rate of mesons is presented for 
proton induced reactions on carbon. 
PACS numbers: 25.75.-q, 25.75.Dw 

1 Introduction 

Currently there is much interest in subthreshold production of heavy mesons 
in relativistic nuclear collisions. It is expected that measurements will re- 
veal the change of the properties of mesons in dense and hot nuclear matter. 
Such predictions are made from theoretical reasons |j] and, indeed, hints to 
a dropping mass of anti-kaons have been found by comparing recent mea- 
surements m with calculations based on transport models ||. For <fi meson 
production there are only a few measurements Q in the reactions of Ni+Ni 
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at 1.93 A-GeV and Ru+Ru at 1.69 A-GeV. Although only a limited amount 
of the phase space was accessible in these experiments, extrapolations to the 
full phase space point to a surprisingly large production cross section. 

First estimates using transport models |j] based on two-particle cross sec- 
tions, derived in a simple meson exchange model ||, seem to underestimate 
the 4> multiplicity. One should however take in to account the lack of precise 
knowledge concerning elementary 7rA, AN and A A collisions. In addition 
multiparticle processes could contribute. Experimentally, such multi-step 
processes can be investigated by accompanying cluster emission in meson 
production [^]. 

In this work we study whether three-body collisions could remarkably 
contribute to the production of <fi mesons in proton-nucleus and heavy-ion 
collisions. At threshold such a mechanism could be dominant because less 
bombarding energy is required if the projectile nucleon interacts with two 
target nucleons instead with a single one. For 'free' nucleons the threshold 
reduces from 2.6 GeV to 1.8 GeV. 

The consideration of production is interesting since also the K~ pro- 
duction proceeds considerably through an intermediate <fi meson. In the 
near-threshold proton-proton collision roughly half of the K~ mesons come 
from the (f) decay H. K~ mesons have been measured by the KaoS collab- 
oration at the heavy-ion synchrotron SIS of GSI Darmstadt in proton-Au 
collisions for a bombarding energy of 3.5 GeV ||, and the ANKE collabo- 
ration envisages the K~ measurements for collisions of protons on various 
light targets at the cooler synchrotron COSY of FZ Jiilich ||10|| . For a firm 
understanding of these reactions a necessary prerequisite is the theoretical 
understanding of the role of the <fi channels. 

We are going to calculate the cross section for the reaction p+2N — > + 
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3N from the invariant amplitude based on tree-level diagrams. In refs. [ |TT| [12 
several diagrams were explored thoroughly for the case of the two-nucleon 
4> production, i.e. N+N — > 0+2N. As shown there, the available data can 
be described by diagrams with internal meson conversion in a 7rp0 vertex. 
Based on this investigation we restrict ourselves to such types of diagrams 
as illustrated in Fig. [p. Three nucleons in states Nx, N 2 , N 3 act together 
via boson exchanges. Literally in the first diagram, after the interaction 
of particles 1 and 2, the intermediate fermion X or pion 7r 2 has gathered 
sufficient energy to produce the meson interacting with the third nucleon. 
In the right hand diagram the time ordering is exchanged. 

Similar processes are described within the current transport models by 
sequential two-step processes which also allow to accumulate the energy of 
several nucleons. In this case, however, the intermediate particles are real, 
i.e. they are moving on-shell. Recently also effort is made to implement 
properly the particles' off-shell propagation (see ref. [fL3]|). In the diagrams 
displayed in Fig. [I] the intermediate particles X and 7r 2 can be either on shell 
or off shell depending on the initial and final momenta of the particles. There 
are subspaces in the phase space, where the intermediate particles move on 
shell. It is our aim to investigate these two effects and to relate them to the 
standard treatment as sequential two-step processes. 

Many-body collisions have already been treated several times in the lit- 
erature. In ref. |TJ] this mechanism was applied to analyse the cumulative 
backward emission of nucleons. In this region of the phase-space off-shell 
contributions become important. Many-body collisions were also incorpo- 



rated in the framework of transport models [|15|, 16]. The results of these 
investigations show moderate effects; here we extend these studies to very 
subthreshold reactions for the production. 
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2 Elementary cross section 



To calculate the cross section for three-particle collisions we will apply the 
one-boson exchange model which is often used to parametrize Lorentz invari- 
ant cross sections. Our approach is strictly based on tree-level diagrams with 
effective parameters adjusted to experimental data. For our investigations 
we restrict ourselves to the diagrams displayed in Fig. [L| We also consider 
two further diagram types, where the p and the 7r 2 mesons are interchanged. 
For the sake of simplicity in the present exploratory investigation we employ 
only the nucleon propagator for the intermediate particle X and the pion 
exchange for the interaction between nucleons. Under these propositions the 
whole set of diagrams we are considering consists of 4 x 36 individual dia- 
grams, where the second factor denotes the number of permutations of the 
respective three fermion lines in the initial and final states. 

In ref. || the interaction Lagrangians for these diagrams have been pre- 
sented which read in standard notation 

C %NN = -^^-^^f^d^Ti , (1) 
£pNN = -gpNNtpYrijjp^ - ^^-ipa^Tipd^py , (2) 
= —ep^^dy *7? . (3) 

These equations contain the nucleon field xp, the iso- vector meson fields if, 
p M and the iso-scalar field The Dirac matrices, the Levi-Civita symbol 
and the corresponding particle masses are denoted by 7 M , e^ap and m^, 
respectively. Greek indices stand for Lorentz indices, and we put h = c = 1. 
Furthermore, a formfactor is needed for each internal line of an exchanged 
meson i coming from or entering a vertex (jkl). We adopt here the standard 
monopole form factor as a function of the square of the transferred four- 
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momentum gf and a cut-off parameter A* H : 



F(2) = ( A )ki? mg 



The following values of the coupling constants and cut-off parame- 
ters have been used: f n NN = 0.989, q p nn = 3.718, f P NN = G-^g P NN, 
= 1.04; A$ NN = 1.6 GeV, A^ = 1.3 GeV, A^ = 0.95 GeV, A£ p0 
= 1.2 GeV. Most of these values are taken from refs. [f|, [lTj and have been 
successfully applied in the near-threshold region. For cf) meson production 
at an energy of 82 MeV above the production threshold we obtain a cross 
section of 0.04 fib using these parameters. This value is smaller than the 
recently measured value ||. However, in ref. [12] it has been demonstrated 
that the agreement between data and calculation is improved by including 
in addition the <fi meson bremsstrahlung from external nucleon legs and the 
final state interaction. 

The cross section for <f) meson production by an incoming nucleon with 
momentum pi = pi^ hitting a nucleon with momentum p 2 which is sur- 
rounded by another nucleon with momentum p 3 in a medium of density n Q 
reads 



AlTTIj^Piq^ 2,TTl]\[ 2 3. (27r) 8 p ro j ec n ons " diagrams 

The first quotient contains the incoming flux with the laboratory momen- 
tum pi ab . The second quotient is related to the density n of the second 
collision partner; the fact that we calculate the cross section at a single nu- 
cleon accounts for the following factor 1/2. The factorial 3! occurs because 
the outgoing nucleons are indistinguishable. The first sum (including 1/8) 
averages over the initial spin-projection quantum-numbers and runs over all 
outgoing spin and isospin states while the second sum contains the 4 x 36 
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diagrams mentioned above which are numerically evaluated including all in- 
terference terms. The integration in Eq. (|^) is carried out over the Lorentz 
invariant phase space spanned by the momenta of the outgoing particles, 
dL mv = U f d 3 p f /(2p° f ) 5 4 (EPf-Ep l ). 



3 On-shell problem 

In some regions of the phase space the intermediate particle X can become 
on-shell and the propagator gets singular. In the following we discuss first 
the situation for the pion propagator ir 2 in Fig. [I] which in a large part of 
the phase space becomes on shell. The standard procedure @, [14], [16], [18|] to 



avoid the singularity is to introduce into the propagator an imaginary part 
imT. That means that one writes the matrix element of Fig. [I] as 



1 

p2 _ m 2 _|_ 



T = f ± 2 ^ I ~^ 2 , (6) 



where the quantities Ti and T 2 describe the two subprocesses 1+2 — > 4+5+7T 
and 7r+3 — > 6+0, which are connected by the propagator of the pion with 
mass m = m n and four- momentum p+ The physical reason for the oc- 
currence of the imaginary part is that the intermediate particle does not 
propagate in vacuum but in a dense medium instead. Then it can move only 
a finite distance of the mean free path A during its life time r which is related 



to the imaginary part in the denominator via T = 1/r (see refs. [H, |16|j). 
The same method is applied if unstable particle are treated in the entrance 
channels I 



However, the on-shell propagation is in fact a process which is already 
treated in standard transport models as two consecutive two-step processes, 
where the particle X is treated as a real on-shell particle. Now we try to 
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separate the genuine three-body processes out of the contributions given by 
the whole propagator in Eq. (|). 

If the particles were in vacuum the amplitude T would also be given 
by Eq. (^), but the quantity mT would be replaced by e with Feynman's 
prescription e — > +0. The square of this matrix element diverges since it is 
proportional to the square of the function 6(p 2 — m 2 ) in the limit e — > 0. 
Since the origin (see e.g. [|IJ) of the S function is the integration over the 
time the pion moves, one can replace it with r5(p 2 — m 2 )/(27rm), where the 
time t is connected with e via r = m/e. Thus, the transition probability 
is proportional to the proper time r in the rest system of the intermediate 
particle. This treatment is analogous to the standard derivation of the 
cross section (|5|) which becomes proportional to the space-time volume. 

Due to the terms coming from the exchange of the initial and final mo- 
menta the matrix T becomes on shell in different regions of the phase space 
spanned by the four final-state momenta. In general there are 3x3 sub- 
spaces defined by (p a — pi,) 2 = m 2 , where p a and p& are the total momenta 
of the three possible initial pairs a =[(12), (23), (31)] and three final pairs b = 
[(4,5), (56), (64)], respectively. 

In the near neighbourhood of each subspace we divide the T matrix in a 
smooth part and a singular one, 



T = T + T 1 - T 2 . (7) 

pZ _ m A _|_ %e 

Contrary to Eq. (|) the amplitudes T\ and T 2 describe now subprocesses with 
free (on-shell) particles. In calculating the phase space integral we treat the 
square of the second term as 

O = CTtwo-step + ®three i (8) 
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7T7~ 

(^two-step = C I dL inv — 5{p 2 - m 2 ) \T\ 2 , (9) 



ni 



/(fp 1 — TTl 2 ^] 2 7T 



(p 2 — m 2 ) 2 7r 
[{p 2 — m 2 

where T = T{r 2 . The quantity C comprises the factors standing in front of 
the integral in Eq. (pj). These expressions arise by splitting the propagator in 
Eq. (0) into the 5 function and its principal value, while the last term a' three 
is constructed such that it is finite as we will see below. 

We identify the first part (|9|) with the two-step process as it is propor- 
tional to the time of flight of the intermediate particle which can only be 
fixed by additional information on the collision geometry. This part can be 
transformed into an expression describing a sequence of two two-nucleon col- 
lisions. Introducing the identity / d A p5 A {p\ + p 2 — Pi — p) = 1 and using the 
definition of the two-particle cross sections we obtain 

O-two-step = J dp n ^~(Pl>P2) 0>(Ptt,P3) {tIqVt) . (11) 

Now the cross section has been factorized in the cross sections da/dp^ 
for pion production and sigma^ for production via pion absorption by the 
nucleon 3. The sum over the final subspaces h has cancelled the factorials 
in the definitions of the cross sections. In principle one should sum Eq. ( [TT] ) 
over the three subspaces a, however, we have used only the first one in 
accordance to Eq. (El). The last factor in the round bracket contains the 



relative velocity v = (pnPzjVreh v re i = y (1 — (mm N / (p n p 3 )) 2 between the 
intermediate particle and its partner 3. It originates from the flux factor 
in the definition of the production cross section and is proportional to the 
length L = vt the pion travels. 

A similar expression like Eq. ( |lTD arises if the intermediate nucleon X 
moves on shell. Formally, the integral over the pion momentum in Eq. (|TT[) is 
replaced by an integral over the scattering angle of the elastic nucleon-nucleon 
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scattering. To derive such an equation the on-shell part of the amplitude T, 



see Eq. @, is transformed into a sum over the spin quantum numbers of the 
intermediate nucleon using the properties of the Dirac operator, 



To arrive at the standard formula the following approximation is needed: 
X) I H n Ti n T2n\ 2 ~ 1/21] \T\\ 2 \T2\ 2 which is fulfilled if spin-flip processes 
are not important. 



number A the cross section would be proportional to A 4 ^ in the subthreshold 
region as the length L is proportional to A 1 ^. However, if the reaction 
proceeds in a dense medium the length reduces to the mean free path, and 
the factor n^vr = \jo to t gives the inverse total cross section. This brings us 
back to the prescription of Eq. (]6|) and the two-step formula employed, e.g., 
in |[20|| . In this case the cross section of a p+A reaction is proportional to A 2 / 3 
taking into account the fact that the incoming particle suffers rescattering 
processes too. 

After separating the two-step process from the total cross section the 
remaining parts form the genuine three-body cross section a three- It contains 
the smooth terms of the T matrix in Eq. (|7|) and the limit of the divergent 
contributions of Eq. (|TUJ) . To see that the expression in the square brackets 
in Eq. (|ToD is finite we analytically integrate over p 2 around the singularity 
within the interval between p 2 = m 2 ± mD. Assuming that T is a smooth 
function of p M we obtain —2/mD for a small but finite value of D. Dividing 
the integration into parts near and far the singularity we write the three-body 
cross section as 





If we apply Eq. ( |TTj ) to a p+A collision of a nucleon on a nucleus of mass 
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-Sip 1 - m 2 ) [ 2mD\T\ 2 + 2nlm(Tf*) - -^-r\f\ 2 

mD 



(13) 



where denotes the step function. Terms of higher than first order in mD 
have been neglected. 

4 In-medium effects 

In Eq. (|T0| ) we have assumed that the intermediate particles move freely. 
Inside a medium we can employ a finite imaginary part as already inserted 
in Eq. Recently great effort has been made to construct the in-medium 
propagators and corresponding spectral functions calculating the self-energies 
in one-loop order and higher approximations to couple the propagating par- 
ticles to the constituents of the surrounding medium. However, for simplicity 
reasons we treat these effects phenomenologically as collision broadening pi[] 
introducing in the propagator in Eq. (H) the collision width T co u = nova to t 

r = - = T (p 2 ) + T coll (14) 
r 

We have already foreseen that the intermediate particle is a resonance with an 
energy dependent decay width r . Although the definitions in Eq. fllO|) can 
be used also for a finite value of e after transforming back the 8 functions 
and their squares into the finite-epsilon representations, this procedure is 
ambiguous here since there is no singularity anymore. Thus, we define the 
genuine three-body cross section as the difference of the total cross sections 
and all possible two-step cross sections: 

Vthree = cj dL mv (\T\ 2 - £ \f c \ 2 5(p 2 c - 7T# ^) . (15) 

The sum runs over all open channels c for two-step collisions, where reso- 
nances are included as long as their decay channels are open, e.g. intermedi- 
ate p mesons with masses m c being larger than the two-pion decay threshold. 
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Due to the surrounding medium the T matrix looses its Lorentz invari- 
ance as the relative velocity refers to the medium's rest frame which we fix 
to be the laboratory system. We mention that the value of v becomes very 
large for pions which create a <ft meson. A pion needs a value of v ~ 11 in 
a collision with a nucleon resting in the laboratory system, while its cross 
section amounts to about 25 mb. The corresponding relative velocity for a p 
meson is considerably smaller due to its larger mass. 

In our calculations we use the effective widths defined in Eq. flH|) for 
those particles only which serve as intermediate particles, and treat particles 
which occur inside the diagrams of the two consecutive processes like free 
particles. This corresponds to the widely used standard method by which 
cross sections are calculated with free particle propagators. 

5 Numerical results 

Now we are going to discuss the cross section for a proton which hits another 
proton at rest embedded in a protonic surrounding with density no=0.16 
fm~ 3 . We vary the bombarding energy Ei^ from the three-body threshold of 
1.8 GeV for <fr meson production up to the two-nucleon threshold of 2.6 GeV. 
In this special case of two particles without relative motion the nucleonic 
intermediate states cannot get on shell below the two-nucleon threshold while 
pions and p mesons already allow two-step processes at an bombarding energy 
of about 0.04 GeV above threshold. 

First we have calculated the cross section without the in-medium effect. 
The two-step cross section (|D|) is proportional to the flight time r or the 
system size. On the other hand the genuine three-body cross section which 
is calculated according to Eq. (|T3| ) does not depend on the system size and 
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is presented in Fig. ^ by the thin full line. 

The cross sections are considerably diminished if the rescattering effects 
of the intermediate particles are taken into account. We have used a total 
cross section of 25 mb for both pions and p mesons. The two corresponding 
two-step processes (defined within the right hand site of Eq. fllBD ) are also 
displayed by the dashed and the dot-dashed lines in Fig. |], respectively. 
These cross sections are considerably smaller than the original three-body 
cross section. The sum of these two quantities has to be compared to the 
total cross section calculated with the complete T matrix of Eq. ( |T5f) which is 
presented by the thick full line. The difference should be the contribution of 
the genuine three-body cross section which, however, turns out to be negative 
(thick dashed line). The reason is that due to the large width of 200 MeV 
the pion propagator in the three-body matrix element reaches out very far 
in the phase space and picks up T matrix elements which are smaller than 
in the on-shell region. 

Finally, we calculate the cross section for the collision of a proton on 
a target nucleus consisting of N neutrons and Z protons. The production 
is now a superposition of primary reactions on pp and pn pairs. The pro- 
duction at a pn pair is strongly enhanced compared with that at a pp pair 
because the isospin coupling allows in the latter case only 7r° exchanges which 
have smaller vertex strengths than those of charged mesons. In more com- 
plete calculations which, e.g., include also the exchange of a mesons such 
a large difference may not occur. Furthermore, the nucleons have Fermi 
momentum which can be taken into account by using the spectral function 
5(p° — m N — p 2 /[2(A — l)m N ] + E sep ) exp (— 5p 2 /2kp ermi ) with a Fermi mo- 
mentum of kFermi = 0.27 GeV and an average separation energy of E sep = 20 
MeV. The effect of the Fermi motion allows the production already at much 
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lower energy. In view of the exploratory nature of our calculations it is not 
worth using improved spectral functions which can be found in ref . [^] . Due 
to the effect of Fermi motion the production already begins at much lower en- 
ergy. In Fig. |3] the cross sections <J V ~ VV and o p ^ pn for <ft production at a proton 
in the surrounding of protons or respectively neutrons are presented. They 
are used to calculate the <fi production in the reaction of a proton with 12 C 
in accordance to a = [(Z(Z — 1)<j p - pp + 2ZNa p - pn + N(N — l)a p - nn ]/A. We 
have used an effective particle number A = 6 as in ref. [[?tj which results from 
the fact that many target nucleons are screened due to scattering processes 
of the incoming proton. For comparison we have also included into Fig. ^ 
the production cross-section (dash-dotted line) calculated as superposition 
of primary two-nucleon collisions, a = Za p _ p + Na p - n , which is smaller than 



the three-body production, in contrast to the results of ref. p0[ . 

Finally we should keep in mind that we have constrained our Lagrangian 
of the nucleon-nucleon-meson interaction. For instance neglecting the in- 
termediate A particles causes an underestimation of the pion production. 
Therefore, the actual <f> production is expected to be larger than presented 
in Figs. H] and [3. 



6 Conclusion 

It has been our aim to study the features and consequences of elementary 
three-body collisions. This is an important issue when calculating particle 
production near threshold for nucleon-nucleus and nucleus-nucleus collisions 
within the frame work of transport models. Here, we have studied the <fi 
production below the free nucleon-nucleon threshold. It turned out that 
genuine three-body processes are not important as most of the intensity can 
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be evaluated via consecutive two-nucleon reactions. In the case considered 
here the two-step processes slightly overestimate the cross section. However, 
it is necessary that all intermediate particles which can become on shell are 
treated properly within the transport codes as they contribute essentially to 
particle production. 
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Figure 1: Dominating tree-level diagrams contributing to the three-body <p me- 
son production. Two further diagrams are used in the text, where the 
7T2 and the p lines are interchanged. 
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Figure 2: Cross section of <fi meson production by a proton hitting a proton at 
rest embedded within a proton environment. The thin line gives the 
three-body cross section without intermediate rescattering whereas the 
other lines show the results in a medium with rescattering. For details 
see text. 
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Figure 3: Cross section for cp meson production for a proton on carbon collision 
(solid curve). Shown are also the cross sections for a proton colliding 
with a proton surrounded by protons (dotted line) or neutrons (dashed 
line) including Fermi motion. The dash-dotted line shows the effect of 
two-nucleon collisions. 
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